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Abstract 

We calculate the energy-band structure of a He atom trapped within the 

interstitial channel between close-packed nanotubes within a bundle and its 

influence on the specific heat of the adsorbed gas. A robust prediction of 

our calculations is that the contribution of the low-density adsorbed gas to 

the specific heat of the nanotube material shows pronounced nonmonotonic 

variations with temperature. These variations are shown to be closely related 

to the band gaps in the adsorbate density of states. 
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1. INTRODUCTION 



The adsorption of gases in and on the bundles of carbon nanotubes has recently attracted 
much attentionS'ii'ilii'i. The carbon nanotube materials provide an excellent opportunity 
for a study of adsorption! in confined geometries. This phenomenon is not only of academic 
interest - a motivation for its study is potential application of the carbon nanotube materials 
for a very efficient storage of gaseJ^'IBI. 

^He atoms adsorbed in strongly bound positions have been detected by measurements of 
desorption from single wall nanotube (SWNT) bundleJl. Recent low temperature specific 
heat measurements on SWNT bundlesili have detected a contribution to the specific heat 
which is very sensitive to the amount of adsorbed gases, in particular ^He. Motivated by 
these experimental results, here we theoretically analyse the behavior of He atoms adsorbed 
in the samples consisting of SWNT bundles. In particular, we investigate the quantum states 
available for the adsorption of He atoms in the interstitial channels between the SWNTs, 
close-packed in a triangular lattice. We analyse the consequences of specificities of energy 
spectrum of adsorbed He atoms on the specific heat pertaining to the adsorbed He gas. 

The effects of interactions between interstitial He adsorbates and their consequences on 
the specific heat have been examined in Ref. ^. Transition of He gas adsorbed in the intersti- 
tial channels to an anisotropic condensed phase has been predicted as a result of attraction 
between atoms in neighboring channels. In this work we concentrate on the interaction of 
He atoms with a surrounding nanotube medium and neglect He-He interactions. 

The main results and the outline of the article are as follows. In Sec. 2, we propose 
a model potential for the interaction of He atoms with surrounding medium. Our model 
of the potential confining He adsorbate to an interstitial channel is sufficiently refined to 
allow the investigation of the infiuence of details of He confinement on the specific heat 
of adsorbate gas. We shall demonstrate that the periodicity of the potential along the 
channel axis causes gaps in the adsorbate density of states. These gaps increase with the 
magnitude of corrugation of the interaction potential along the interstitial channel axis. As 
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a consequence, in Sec. 3 we shall show that the specific heat pertaining to the adsorbate gas 
shows nonmonotonic behavior with temperature. The precise dependence of specific heat 
on the temperature will be shown to be determined by the corrugation of the interaction 
potential. 

Our approach is similar in spirit to the one presented in Ref. ITB|, where the specific heat 



-y^ + V{x,y,z) 



"^^{x^y^z) = E^^^{x,y,z), (1) 



of low-density He films on graphite has been examined. The most important results of this 
article are summarized in Sec. 4. 

2. ADSORPTION POTENTIAL AND BAND STRUCTURE CALCULATION 

The Schrodinger equation for the wavefunction of a single adsorbate within 

the interstitial channel is 

A' 

"2M 

where M is the mass of adsorbate and k denotes a set of quantum numbers needed for the 
specification of the quantum state of the adsorbate. 

To find the quantum states supported by equation (p one needs to construct the potential 
V{x,y,z). However, the construction of the relevant potential is somewhat problematic for 
two reasons. First, the potential describing the interaction of He atom with the carbon 
nanotube in not known with great precision and often the model potentials based on the 
assumption of pairwise additivity of binary He-C interactions are used&i. Second, the details 
of arrangement of nanotubes within a bundle are also not known. In particular, the relative 
axial alignment of nanotubes (relative axial offsets) within a bundle is not known! . As shown 
in Ref. 0, different nanotube alignments can produce very different potentials for adsorption 
of He atom in the interstitial channels. It is also likely that bundles contain a mixture of 
tubes with different wrapping angles0lli. Details of this mixture may vary depending on 
a particular bundle in the sample. Furthermore, the distribution of nanotube diameters 
in the sample has a finite widthlii0. All these facts imply that a detailed information on 
the potential is at present not available and that the relevant potential depends on the 
interstitial channel in question. 



Here we propose a model potential V{x,y,z) which captures all the essential physics 
involved, yet simple enough to allow an efficient and transparent calculation of the quantum 
states and predict the dependence of specific heat on temperature. According to references 
^ and ^, the effective potential confining the motion of He atom in the plane perpendicular 
to the channel axis can be represented by a paraboloid-like shape with the minimum at 
the geometrical center of the channel (see Fig. 1 of Ref. §). The potential varies along the 
channel axis, reflecting the variations of the nanotube electronic density in this direction!. 
The simplest potential which displays the described behavior is 

V{x, y, z) = D + a{x^ + y"^) + 2Vg{x, y) cos{gz) 

= Vo{x, y) + 2Vg{x, y) cos{gz), (2) 

where g = 27T/a is the inverse lattice vector associated with the periodicity of the potential 
along the channel axis {z direction) and the factor of 2 in front of Vg{x,y) term has been 
introduced for convenience. The origin of x and y axes is positioned at the center of the 
channel and the channel axis is perpendicular to the xy plane. D defines the minimum 
of the potential at the center of the channel, in absence of the corrugation term, Vg{x,y). 
Equation (|T]) can be solved by writing the wavefunction as 



*^(x,y,z)= I^m,n,k+G'^rn,n,k+G{x,y,z), (3) 

m,n,G 



where 



<^^,n,k+G{x, y, z) = f HUSx)H^{5y) exp (-ls\x' + y')) e^^'+''> 

= \m,n;k + G > (4) 

is the solution of the hamiltonian containing only Vo{x,y) part of the potential. Here, is 
the box quantization length along the axis of the channel, Hm,n are the Hermite polinomials, 
(5^ = 2Ma/h'^, m and n are the harmonic oscillator quantum numbers, k is the wavevector 
associated with the motion of He atoms in z direction, and G is a reciprocal lattice vector 
given as G = pg., with p integer, associated with the periodicity of the potential in z direction. 
The Schrodinger equation (0) reduces to matrix equation 



E 

m',n\G' 



+ <m,n;k + G\2Vg{x,y)cos{gz)\m',n';k + G' > (3^^^^^,^^G = ^- (5) 

m',n',G' 

Em',n' are the eigen energies of the hamiltonian without the corrugation part of the potential 
and are given as 



Em',n' = D + hujo{m +n + 1), 



(6) 



where luq = y 2a/M is the characteristic frequency of the harmonic oscillator. The spectrum 
of energies of adsorbed He atom can be found by looking for nontrivial solutions of equation 
(^). For fixed k, this produces a number of different solutions (bands) denoted by index A. 

By systematic examination of a range of k values, it is possible to obtain the full band 
structure, Ex{k). Similar approach has been used for a calculation of band structure of He 
atoms adsorbed on the surface of graphiteEl and NaCl@. 

We parametrise the corrugation term in the potential as 



Vg{x,y) = co + C2(a;^ + ?/^ 



(7) 



The matrix elements of the corrugation term can be cast in the analytic form as 



V^,^Q^ = < m,n; k + G\2Vg{x, y) co?,{gz)\'m' , n'; k + G' > 



+ 



_C2_ 



Co + ^{m' + n' + 1] 



52 



5n,n' {m' + l)(m' + 2)5m,m'+2 + \J m' {m' - l)5m,m'-2 



(^^J{n' + l){n' + 2)5n,n'+2 + ^n'{n' - l)5n,n'-^'\ } • (8) 

The parameters of the proposed model of interaction potential can be fixed by examining the 
potential obtained as a pairwise sum of binary He-C site interactions. We model the effective 
He-C site interaction with a Lennard- Jones (LJ) form, 4e[(cr/r)-'^2 _ (o-/r)^], where r = |r| 
denotes the distance between He atom and a C site. Various different forms of effective 
He-C potential, based on the analysis of He-graphite interaction, have been proposed in 
the literature0'0. The most important deficiency of the LJ form is the assumption that 



the potential depends only on the magnitude of r0. However, in the view of a lack of 
knowledge on the precise geometry of a bundle, the LJ form should serve well in obtaining 
the main features of the interaction potential. The LJ parameters we adopt are e = 1.34 



meV, a = 2.75 A. This choice of parameters was suggested in Ref. |T9|, where the LJ form of 
the potential has been optimized to reproduce the experimental data on ^He scattering from 
the basal (1000) plane of graphite. These parameters are slightly different from the ones 
reported in Ref. |^ (e = 1.46 meV, a = 2.98 A) which were obtained from the semiempirical 
combining ruleJ^. 

In the following, we consider (10,10), "armchair" carbon nanotubes0 with a diameter of 
d = 13.8 A. The centers of the nanotubes are separated by 13.8 A+3.2 A= 17.0 A, where the 
value of the intertube separation of 3.2 A has been adopted, in agreement with experimental 
findings^ and also with Refs. ^ and ^ Assuming that all the tubes surrounding a channel 
are axially aligned, and that the whole structure of the bundle can be obtained from a single 
tube by applying a two dimensional translation characterised by a general two dimensional 
lattice vector given as t = iti + jt2, where i and j are integers and ti and t2 are basis 
vectors of a triangular lattice of the bundle, the parameters consistent with our model of 
the potential are D = -36.8 meV, a = 44.7 meV/A^, cq = -0.05 meV, C2 = 10.7 meV/A^ 
and g = 27r/2.5 A= 2.51 1/A. This set of parameters provides an excellent fit to the angular 
average of the total potential^. 

The density of states per unit length of the channel can be obtained from 

9{E) = ^j:6{E-E,{k)). (9) 

X,k 

Some care should be taken in evaluation of Eq. (|) to properly account for the degeneracy 
present in the energy spectrum, E\{k). In Fig. |1| we present the band structure and the 
corresponding density of states obtained from the determined set of potential parameters. 
For convenience of presentation, the delta functions in equation (^) have been broadened 
to gaussians of width 0.03 meV. Note the enhancement of the lowest band effective mass of 
He atom, M* = 2.37M = 9.5 amu and the band gap between -27.75 meV and -25.5 meV. 
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The bandwidth of the lowest energy band is 0.19 meV (2.2 K). This should be compared 
with the value of 0.18 K found for (18,0) tubes in Ref. |^. The difference between these two 
values is determined in part by the different periodicities (cell lengths) of the "zigzag" and 
"armchair" carbon nanotubes. 

The mass enhancement should be compared with the one found by authors of Ref. 
M* = 1.3M. The difference between their and our value can be explained by somewhat 
larger intertube separation presumed in the calculations in Ref. || (3.3 A, whereas we used 
3.2 A). As shown in Ref. ^this parameter strongly influences the magnitude of the potential 
corrugation and thus also the effective mass enhancement. The ground-state energy, Er, 
of adsorbed *^He atoms as found from our model potential is -27.94 meV. The ground-state 
energy found in references ^ and |^ for nanotubes of the same diameter (13.8 A), but for the 
interaction potential which does not exhibit corrugation along the channel axis was -29.07 
meV and -33.3 meV, respectively. This shows consistence of our potential with the potentials 
previously proposed in the literature. The experimental estimate of the binding energy based 
on desorption measurementsi was reported to be -28.5 meV, in good agreement with the 
value we obtain. The subsequent reevaluation0 of the experimental data first presented in 
Ref. 1^ yielded a smaller binding energy (-19.8 meV). The authors of Ref. suggested that 



this value of binding energy may be compatible with adsorption of He on the outer grooves 
of the bundle. Our results should be of help to determine the relevant adsorption sites. 

One of the attractive features of our potential model is that the various different arrange- 
ments of nanotubes surrounding a channel can be easily (but approximately) simulated by 
changing the parameters of the potential, especially those representing the corrugation of 
the potential. Additionally, the effects of possible inadequacies of the binary He-C site po- 
tential (such as the lack of anisotropy in the effective interaction) can be avoided by treating 
the parameters Cq and C2 of the potential as variables. In Fig. ^ we present the behavior 
of the density of states as a function of parameter C2 which defines the magnitude of the 
potential corrugation. As expected, for smoother potentials, the band gaps, binding energy 
and effective mass are smaller. In particular, the lowest band mass enhancements are 2.37, 



1.23 and 1.02 for the calculations presented in panels a), b) and c) of Fig. respectively. 
3. SPECIFIC HEAT OF LOW DENSITY ^He INTERSTITIAL GAS 



Knowing the density of states, we can calculate the isosteric specific heat of He atoms 
adsorbed in the nanotube interstitial channels. The approach we pursue is similar to the 
calculation of specific heat of the adsorbed (sub)monolayer of noninteracting atoms on a 
corrugated surface discussed in Ref. |1^. The influence of adsorbate-adsorbate interactions 
on the specific heat of the He overlayer on graphite has been studied in Ref. 

The total energy of a system of N noninteracting ^He atoms is given by 

roc 

U = lJ g{E)Ef{E,T)dE. (10) 
The Bose-Einstein distribution function in the above equation is given by 

mT)= \ (11) 

where fc^ is the Boltzmann constant and T is temperature. The chemical potential /i can 
be determined from the requirement of the conservation of number of adsorbates, N , 



N = LJ g{E)f{E,T)dE. (12) 

JEr 

Solution of Eq. (|12D results in the temperature dependence of chemical potential /i for a 
fixed total number of ^He atoms. The isosteric specific heat C/N can be obtained from 



where integral quantities Ij are defined as 

'E-^i{T) 



I,= / g{E)E^exp 



fiE,T)dE, J = 0,1,2, (14) 



in accordance with the notation of Ref. |1^. Equations (|T3|) and ( p!^ are very convenient for 
numerical implementation once the dependence of chemical potential on temperature, fi{T) 
has been found from Eq. i\L^ . 
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In Fig. ^ we plot the isosteric specific heat as a function of temperature for three 
different average adsorbate linear densities and with potential parameters obtained from the 
assumption of perfect alignment of the tubes. The dashed (full) line corresponds to average 
interadsorbate distance of 1 nm (2 nm). It is interesting to note here that the mean He-He 
distance along the channel estimated from experimental data in Ref. |^ was found to be 
about 3 nm. The dash-dotted line in Fig. |^ corresponds to this average linear ^He density. 

As noted in Ref. |], perfectly aligned tubes result in the highest corrugation of the po- 
tential seen by He adsorbate. In this respect, it is interesting to see how the specific heat 
changes as a function of a potential corrugation along the channel axis. In Fig. ^ we plot 
the variation of specific heat with temperature for three different corrugation parameters 
C2 and we fix the average linear density of He adsorbates to N/L^ = 0.1 1/A. All three 
curves display a nonmonotonic behavior and a minimum in specific heat around T = 4 K 
is observed. At higher temperatures all curves approach to the value of 0.5 (thin dotted 
line), refiecting a specific heat characteristic of a particle with one translational degree of 
freedom. For even higher temperatures (~ 20 K, not shown in the figure) the excitation of 
higher harmonic oscillator levels [m! or n' in Eq. (^) different from zero] becomes possible 
and the specific heat becomes larger than 0.5. For such high temperatures, desorption of 
He from the sample becomes probable0'i. The minimum in specific heat occurs due to the 
existence of the band gap in the density of statesElil. 

As the band gap becomes smaller (i.e. for lower corrugations of the potential) the min- 
imum in specific heat becomes less pronounced and practically disappears for the smallest 
corrugation of the potential considered in the calculations presented in Fig. H. It is inter- 



esting to note that our results show similar trends as the ones presented in Ref. |T^. Due 
to the different dimensionality in our case, the specific heat is a factor of about 2 smaller 
than the one calculated for rare overlayer of He atoms on graphitei. Also, due to the ef- 
fective magnification of the corrugation in a restricted geometry of the interstitial channel, 
the variation of the specific heat with temperature is significantly more pronounced. 
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4. SUMMARY AND CONCLUSIONS 



A robust prediction of our model is appearance of band gaps in the density of states of 
He adsorbates in the interstitial channels. This prediction is based on the assumption of 
periodicity of the potential along the channel axis and it holds irrespective of the precise 
alignment of the tubes surrounding a channel, as long as all three tubes have the same 
internal symmetry. The existence of the band gap causes a nonmonotonic behavior of 
the specific heat associated with the interstitial He gas, an effect previously observed for 
low-density overlayer of He atoms on a graphite surfaceEl. We observe an appearance of 
the minimum in specific heat around T = 4 K. As the corrugation of the potential along 
the channel axis smoothens, the minimum in the specific heat becomes less dramatic and 
ultimately vanishes. The effects we predict should be easily tested experimentally in the 
regimes where the specific heat of the adsorbates can be easily distinguished from the overall 
specific heat (for low temperatures and not too low concentrations of He adsorbatesIS) . 
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FIGURES 

FIG. 1. a) Calculated band structure of a ^He atom adsorbed in the interstitial channel, b) 
The corresponding density of states per unit length. 

FIG. 2. Density of states of "^He atoms adsorbed in the interstitial channel as a function of the 
magnitude of corrugation of interaction potential (c2 parameter), a) cq = —0.05 meV, C2=10.7 
meV/A^. b) Co = -0.05 meV, C2=5 meV/A^. c) cq = -0.05 meV, C2=2 meV/A^. 

FIG. 3. Isosteric specific heat of ^He gas within the nanotube interstitial channels as a function 
of the sample temperature for three different average He linear densities, N/Lz- Dash-dotted line: 
A^/L^=0.033 1/A. Full line: 7V/L^=0.05 1/A. Dashed line: Ar/L^=0.1 1/A. 

FIG. 4. Dependence of specific heat of ^ He gas on the corrugation of the interaction potential 
along the channel axis. Full line: cq = —0.05 meV, C2=10.7 meV/A^. Dashed line: cq = —0.05 
meV, C2=5 meV/A^. Dash-dotted line: cq = —0.05 meV, C2=2 meV/A^. In these calculations, 
average linear density of ^He atoms was set to 0.1 1/A. 
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